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Abstract 

Background: The green crab Carcinus maenas is known for its high acclimation potential to varying environmental 
abiotic conditions. A high ability for ion and acid-base regulation is mainly based on an efficient regulation 
apparatus located in gill epithelia. However, at present it is neither known which ion transport proteins play a key 
role in the acid-base compensation response nor how gill epithelia respond to elevated seawater pC0 2 as 
predicted for the future. In order to promote our understanding of the responses of green crab acid-base 
regulatory epithelia to high pC0 2 , Baltic Sea green crabs were exposed to a pC0 2 of 400 Pa. Gills were screened 
for differentially expressed gene transcripts using a 4,462-feature microarray and quantitative real-time PCR. 

Results: Crabs responded mainly through fine scale adjustment of gene expression to elevated pC0 2 . However, 2% 
of all investigated transcripts were significantly regulated 1.3 to 2.2-fold upon one-week exposure to C0 2 stress. 
Most of the genes known to code for proteins involved in osmo- and acid-base regulation, as well as cellular stress 
response, were were not impacted by elevated pC0 2 . However, after one week of exposure, significant changes 
were detected in a calcium-activated chloride channel, a hyperpolarization activated nucleotide-gated potassium 
channel, a tetraspanin, and an integrin. Furthermore, a putative syntaxin-binding protein, a protein of the 
transmembrane 9 superfamily, and a CI7HC0 3 " exchanger of the SLC 4 family were differentially regulated. These 
genes were also affected in a previously published hypoosmotic acclimation response study. 

Conclusions: The moderate, but specific response of C. maenas gill gene expression indicates that (1) seawater 
acidification does not act as a strong stressor on the cellular level in gill epithelia; (2) the response to hypercapnia 
is to some degree comparable to a hypoosmotic acclimation response; (3) the specialization of each of the 
posterior gill arches might go beyond what has been demonstrated up to date; and (4) a re-configuration of gill 
epithelia might occur in response to hypercapnia. 



Background 

With increasing atmospheric pC0 2 , a decrease in global 
surface ocean pH of between 0.4 to 0.8 units is pre- 
dicted due to oceanic CO2 uptake [1,2]. Although large 
changes in atmospheric C0 2 have been recorded 
throughout earth history, the current anthropogenic 
increase in pC0 2 is much more rapid and severe than 
the cyclic changes of pC0 2 during the last 20 million 
years [3]. The resulting changes in carbonate chemistry 
speciation, termed 'ocean acidification', may become a 
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general stress factor modulating future marine commu- 
nities by differentially influencing the fitness of marine 
species [4-7]. Many studies suggest complications for 
marine ectothermic metazoans in response to future 
ocean acidification: reduced growth and calcification 
rates, reduced rates of development, altered energy bud- 
gets, disturbed acid-base status, disturbed chemosensory 
function and even increased rates of mortality have been 
measured [4,7-14]. On the other hand, increases in 
growth and calcification rate were shown primarily in 
organisms that accumulate significant concentrations of 
bicarbonate in their body fluids, e.g. decapod crusta- 
ceans [12], cephalopods [15] and fish [16]. 



O© 2011 Fehsenfeld et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BiolVlGCl C^ntrBl Commons Att ibution License (http://creativecommons.0rg/licenses/by/2.O), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Fehsenfeld et al. BMC Genomics 201 1, 12:488 
http://www.biomedcentral.eom/1 471 -2 1 64/1 2/488 



Page 2 of 1 7 



Coping with ocean acidification 

Biological impacts of acidification are strongly related to 
life history, genetic pre-disposition and physiological 
acclimation potential of the species in question. Elevated 
environmental pC0 2 results in an increased extracellular 
C0 2 partial pressure, as positive diffusion gradients of 
C0 2 have to be maintained in order to excrete meta- 
bolic CO2 [5]. This can then lead to an acidification of 
extracellular fluids [17-19]. However, several active, high 
metabolic species with pH sensitive respiratory pigments 
regulate extracellular pH (pH e ): active modulation of the 
extracellular carbonate system leads to bicarbonate 
accumulation and pH compensation while maintaining 
pC0 2 values sufficiently high for diffusive C0 2 flux out 
of the animal [fish: 20; crustaceans: 21; cephalopods: 22, 
23]. In teleost fish, cephalopods, and decapod crusta- 
ceans, the majority of the acid-base relevant ion regula- 
tory apparatus is located in gill epithelia. It is thought 
that net proton extrusion is primarily achieved via active 
(V-type H + -ATPase) and secondarily active ion transport 
molecules (e.g. sodium proton exchangers, NHE; sodium 
bicarbonate cotransporters, NBC), with a strong sup- 
porting role of carbonic anhydrases (CAs) and Na + /K 
+ -ATPase [24]. 

Carcinus maenas 

In the case of the hyperosmoregulator Carcinus maenas 
(green or shore crab), the three posterior gills number 
7-9 have been found to be involved in osmo- and acid- 
base regulation, while the anterior six gill pairs serve the 
demands of gas exchange [25-27]. Studies of Truchot 
[25] demonstrated that C. maenas actively and rapidly 
accumulates HC0 3 ~ in its hemolymph during acute 
exposure to elevated seawater pC0 2 (690 Pa). However, 
active adjustment of the extracellular carbonate system 
might not only be a short-term response in this species, 
but a persistent physiological response in order to main- 
tain high extracellular pH (pHe). Despite efficient pHe 
regulatory control, other physiological aspects might be 
negatively impacted by hypercapnia and might show 
effects on the animal's performance only in the long 
run. Up to date, only few studies investigated long-term 
effects of acid-base disturbance in crustaceans on the 
physiological level. Potential reponses may include an 
altered calcification rate, as demonstrated for three crus- 
tacean species (the crab Callinectes sapidus, the shrimp 
Penaeus plebejus and lobster Homarus americanus) by 
Ries et al. [12], and other invertebrates, such as the cut- 
tlefish Sepia officinalis [15]. Other adaptations to long- 
term hypercapnia may be similar to what has been 
observed as responses upon hypoxia, both factors being 
closely linked to each other. These might include energy 
conserving strategies (e.g. down-regulation of protein 
synthesis and modification of certain regulatory 



enzymes) as can be observed in response to hypoxia in 
mussels [28], fish [29] or reptiles [30] (reviewed by [31]). 

On the molecular level, it is not known at present 
which ion transporters play a key role in the short- and 
the long-term response to hypercapnia in decapod crus- 
taceans. Although models for the organization of the gill 
epithelia of euryhaline crabs have been postulated 
[32,33], the transporter inventory in decapod crustacean 
gill epithelia and their functional interactions are not 
fully understood at present. As C. maenas is exposed to 
regular short-term fluctuations in pH in its highly vari- 
able habitat in the western Baltic Sea [19], not only 
changes in gene expression levels may play a role in its 
acid-base compensation response. A compensatory 
stress response can as well take place on the post-tran- 
scriptional level e.g. through covalent modification of 
certain regulatory enzymes and their activity through 
phosphorylation-dephosphorylation reactions. This has 
been shown e.g. for the enzyme phosphofructokinase, 
the key enzyme in controlling anaerobic carbon flow in 
glycolysis in mussels in response to hypoxia [34], and 
fish [35]. A post-translational alteration has also been 
hypothetized for the enzyme carbonic anhydrase in the 
blue crab Callinectes sapidus, a close relative of C. mae- 
nas [36]. 

Aim of the study 

In order to promote our understanding of ion and acid- 
base regulatory processes, green crabs from the Baltic 
Sea were exposed to control (39 Pa) and elevated (400 
Pa) pC0 2 (100 Pa = 987 uatm). The experimental condi- 
tions in the laboratory were chosen to mimick the para- 
meters in the crabs' natural environment, the brackish 
water system of the western Baltic Sea. The chosen 
experimental level of pC0 2 enrichment (400 Pa) repre- 
sents a stress scenario that could occur in the natural 
habitat within the next 100 years: pH in Kiel Fjord can 
reach values as low as 7.5 (corresponding pC0 2 = 230 
Pa) in surface waters during summer and autumn. 
Future changes in pC0 2 will be more severe in this 
habitat than in the average surface ocean [19]. As a 
strong extracellular acid-base reaction has been 
observed in C. maenas exposed to hypercapnia [25], we 
hypothesize that the ion regulatory transcriptome of C. 
maenas will respond to hypercapnia with a short- and 
long-term adjustment of expression of important ion- 
and acid-base transporter candidates (e.g. Na + /K 
+ -ATPase, Na + /H + exchanger (NHE), V - type H 
+ -ATPases, CT/HCO3' exchangers and carbonic anhy- 
drases; [33]). As little is known about candidate genes 
for acid-base regulation in crustaceans, we chose to uti- 
lize a microarray approach to screen for previously 
unrecognized candidates. Therefore, gene expression 
profiles of gills of short-term (3 and 7 days) exposed 
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crabs were investigated using a 4,462-feature microarray 
assay recently developed by Towle et al. [37]. Expression 
levels of distinct candidate genes were further investi- 
gated using quantitative real-time PCR analysis following 
long-term exposure to hypercapnia (11 weeks). Expres- 
sion profiles obtained in our study were compared to 
those obtained in a recent hyposmotic study, in which 
gene expression changes after transfer of C. maenas 
from a salinity of 32 to 15 was investigated using the 
same microarray [37]. A second focus of interest was 
placed on genes known to be in involved in the cellular 
stress-response (e.g. heat-shock proteins, superoxide dis- 
mutase, glutathione reductase [38]). 

This is the first study that elucidates the impact ele- 
vated seawater pC0 2 has on the gill transcriptome of a 
decapod crustacean. 

Methods 

Animals, exposure, and tissue sampling 
Short-term C0 2 experiments 

Green crabs (Carcinus maenas) were caught in traps in 
April 2009 in Kiel Fjord (Baltic Sea) at the IFM-GEO- 
MAR pier (westshore building) at 3-4 m depth (54° 
19.8'N; 10°9.0'E). In order to acclimate the animals to 
laboratory conditions, they were kept in a flow-through 
tank (H30 x W30 x L50 cm) in a climate chamber at 
IFM-GEOMAR with ambient aerated Baltic Sea brackish 
water at 13°C for 3 to 7 days. Experimental animals 
were fed ad libitum with mussels (Mytilus edulis) dur- 
ing this time. Female crabs (n = 24, 6 animals per CO2 
treatment and day) with a mean carapace width of 4.6 ± 
0.2 cm were then chosen for the short-term experiment 
and two each were transferred to 20 L tanks of a flow- 
through seawater C0 2 manipulation system (flow rate = 
100 mL min' ) that was supplied with ambient Baltic 
Sea brackish water of a salinity of 14 - 15 (Practical sali- 
nity scale, as it will be used throughout the script). 
Temperature was held constant in the storage tank by a 
thermostat set to 13°C. A light-dark cycle of 12:12 h 
was established. The system design was identical to that 
presented in Thomsen et al. [19]. 

After 24 hours of acclimation in the experimental sys- 
tem resembling the natural parameters encountered in 
Kiel Bight (S = 14.9, T = 12.9 ± 0.1°C), crabs were 
exposed to control (53 Pa, pH = 8.11 ± 0.04) and ele- 
vated (440 Pa, pH = 7.24 ± 0.03) pC0 2 for 3 or 7 days. 
C0 2 was provided by a central automatic C0 2 -mixing- 
facility (Linde Gas, HTK Hamburg, Germany). Animals 
were fed ad libitum with crushed mussels {Mytilus 
edulis). 

After 3 days, one of the two animals in each tank was 
removed for analysis, anaesthetized on ice and then 
killed by destroying its ventral ganglion and removal of 
the carapace. Left gills #7 and #9 were carefully removed 



with forceps and immediately stored in RNAlater" 
(Ambion, #AM7021) until microarray analysis at the 
Mount Desert Island Laboratory, Salisbury Cove, Maine/ 
USA, or quantitative real-time PCR (qRT-PCR at the 
IFM-GEOMAR, Kiel. The procedure was repeated for 
the second animal of each tank after 7 days of 
incubation. 

For qRT-PCR on gills of short-term exposed animals, 
a repeated incubation of Carcinus maenas specimens 
was performed in May 2010. Again, Baltic Sea green 
crabs were caught in traps in Kiel Fjord at the IFM- 
GEOMAR pier. Females had a mean carapace width of 
4.8 ± 0.6 cm. System parameters in the acclimation 
phase were S = 13.4 and T = 10.8 ± 0.1°C. The animals 
(n = 12, 6 animals per C0 2 treatment and day) were 
then exposed to seawater of pH 7.99 ± 0.02 (68 Pa C0 2 ) 
and pH 7.28 ± 0.05 (361 Pa C0 2 ) as described above. 
Gills were isolated only after 7 days in this experiment. 
Abiotic seawater parameters 

To evaluate the carbonate system parameters, water 
samples (500 mL) were taken of 4 randomly chosen 
tanks each for day 3 and day 7 of the short-term hyper- 
capnia experiment in 2009. For the second short-term 
exposure in 2010, water samples of 4 control and 4 
experimental aquaria were taken only on day 7. Total 
dissolved inorganic carbon (C T = [C0 2 *] + [HC0 3 ~] + 
[CO3 2 ]) was measured according to Dickson and Mill- 
ero [39]. In case of the second 1 week exposure in 2010, 
Ct was measured using an AIRICA CT analyzer (Mar- 
ianda GmbH, Kiel, Germany). In all experiments, sali- 
nity, pH and temperature were assessed daily using a 
pH- and salinometer (WTW 340i pH-analyzer/WTW 
SenTix 81-measuring chain, WTW cond 315i/WTW 
TETRACON 325-measuring chain). 

Seawater total alkalinity C4 T ) and pC0 2 was calculated 
from pH and C T using C02SYS software [40] and the 
appropriate parameter and constants (dissociation con- 
stants Kl and K2 according to Roy et al. [41], KHSO4 
dissociation constant after Dickson [42], NBS scale 
[mol/kg H 2 0]; table 1). 
Long-term C0 2 experiments 

Green crabs {Carcinus maenas) for the long-term 
experiment were part of a study by Appelhans et al. 
(under review) that was conducted in parallel. Experi- 
mental animals for this long-term study were caught in 
traps in April 2009 in Kiel Fjord (Baltic Sea) at the IFM- 
GEOMAR pier (westshore building) in 3 to 4 m water 
depth (54°19.8'N; 10°9.0'E). Animals with a carapace 
width of 4.5 ± 0.5 cm were chosen irrespectively of their 
sex (n = 9 animals per C0 2 treatment). Crabs were 
exposed to control (40 Pa, pH = 8.06 ± 0.01) and ele- 
vated pC0 2 (230 Pa, pH = 7.36 ± 0.01). Animals were 
fed with 2 mussels (size = 2.50 ± 0.35 cm) three times a 
week. To evaluate the carbonate system parameters, 
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Table 1 Carbonate system parameters of all hypercapnia experiments with Carcinus maenas 





1 week experiment 2009 


1 week experiment 2010 


1 1 week experiment 2009 


Parameter 


control 


experiment 


control 


experiment 


control 


experiment 


C T [umol/kg] (measured) 


2025 ± 29 


2242 ± 37 


1879 ± 7 


2047 ± 1 0 


1 999 ± 1 8 


2209 ± 24 


A f [umol/kg] (calculated) 


2098 ± 29 


2062 ± 26 


1904 ± 5 


1887 ± 9 


2046 ± 1 9 


2061 ± 20 


pH (measured) 


8.12 ± 0.02 


7.24 ± 0.02 


8.00 ± 0.02 


7.28 + 0.01 


8.06 ± 0.01 


7.36 ± 0.01 


S (measured) 


15.35 ± 0.06 


15.30 ± 0.00 


13.40 + 0.02 


1 3.4 ± 0.02 


14.80 ± 0.13 


14.80 ± 0.13 


T (measured) [°C] 


12.90 ± 0.16 


1 2.98 + 0.05 


10.87 ± 0.02 


10.77 ± 0.13 


1 2.90 + 0.07 


1 2.90 + 0.07 


pC0 2 (calculated) [Pa] 


53 ± 2 


440 ± 25 


68 ± 3 


361 ± 9 


40 ± 5 


230 ± 21 



Measured parameters of the carbonate system were assessed daily during the acclimation phase of 1 or 1 1 weeks, while A T and pC0 2 were calculated applying 
the measured parameters and C02SYS software (see Material and Methods for details). Under experimental conditions (elevated pC0 2 ), a drop in pH of 0.7 - 0.9 
units could be observed. The carbonate system remained stable during the acclimation phase in both, short- (1 week) and long-term (1 1 weeks) experiments. C T 
= total dissolved inorganic carbon, A T = total alkalinity, S = salinity, T = temperature, pC0 2 = partial pressure of C0 2 . Values are given as mean with standard 
deviation. 



water samples of three aquaria per treatment level were 
taken at the beginning, the intermediate phase and at 
the end of the experiment, respectively. Analysis of C T 
was coulometric (SOMMA System autoanalyser, Mar- 
ianda GmbH, Kiel, Germany) and A T analyzed via 
potentiometric titration (VINDTA autoanalyser, Mar- 
ianda; table 1). An overview over the different experi- 
ments and sampling points is depicted in Figure 1. 

Microarray experiment 
RNA extraction 

RNA from the first short term experiment (2009) was 
extracted from whole left gills #7 and #9 (control PCO2, 



Short-term exposure of Carcinus maenas Long-term exposure of Carcinus maenas 
to control (39 Pa) + elevated (400 Pa) pC0 2 to control (39 Pa) + elevated (400 Pa) pC0 2 

(2009) 



(2009) (2010) 

i I 

Gill sampling Gill sampling 

/\ 

day 3 \ 
(gill 7 + 9) \ 

day 7 day 7 

(gill 7 + 9) (gill 9) 



RNA isolation RNA isolation 



Gill sampling 



11 weeks 
(gill 7) 



±± 

Microarray analysis 



qRT-PCR 



Figure 1 Overview over the experiments performed on 
Carcinus maenas specimen of the Baltic Sea A short-term 
incubation experiment was conducted in April 2009 and microarray 
analysis was performed on isolated RNA of gills of these animals 
exposed to control (nominal 39 Pa) and elevated (nominal 400 Pa) 
pC0 2 for 3 and 7 days. The experiment was repeated in 2010 and 
this time, RNA was isolated for quantitiative real-time polymerase 
chain reaction (qRT-PCR) of gill 9 after 7 days. In parallel, a long- 
term experiment was performed on Baltic Sea C. maenas in 2009 by 
Appelhans et al. (under review), and RNA of gill 7 was isolated and 
analysed by qRT-PCR. 



n = 6; elevated pC0 2 , n = 6) using the RNeasy Midi-Kit 
from Qiagen (#75144). Gills were homogenized with an 
OMNI international TH/G7-195STW for 50 s in 4 mL 
buffer RLT from the Qiagen kit. Extraction then fol- 
lowed the manufacturer's description. RNA quality was 
monitored by electrophoresis via the Agilent 2100 Bioa- 
nalyzer system. Each RNA sample was then reversely 
transcribed into cDNA and in parallel fluorescence- 
labeled using the Superscript™ Plus Direct cDNA 
Labeling System from Invitrogen (#L1015-06), following 
the provided protocol. Using both anchored oligo(dT)20 
and random hexamer primers included in the kit, con- 
trol cDNA was labeled with AlexaFluor555 (green fluor- 
escent, equivalent to Cy™3; excitation: 555 nm, 
emission: 565 nm) while cDNA of the C0 2 -treatment 
samples was labeled with AlexaFluor647 (red fluores- 
cent, equivalent to Cy™5; excitation: 650 nm, emission: 
670 nm). 

Microarray assay for Carcinus maenas 

The C. maenas microarray was developed by D. W. 
Towle et al. and prepared as described in Towle et al. 
[37]. The annotation of the respective contigs included 
in the microarray, based on a Carcinus maneas EST 
library (deposited at http://www.ncbi.nlm.nih.gov), was 
performed by Towle et al. [37] using BLASTX analysis 
on a local TimeLogic DeCypher server (Active Motif, 
Inc.). The most informative BLASTX hit for each contig 
sequence was selected manually from the ten highest 
scoring hits identified by the programm. In order to 
improve the information on each included gene, a re- 
annotation of the contigs has been performed in 2009 
by the authors of this study and compared to the exist- 
ing annotation (Additional File 1 Table SI). A sequence 
blast with the online portal Blast2GO [43] was per- 
formed to improve the quality of the annotations. Addi- 
tionally, a mapping and substance Gene Ontology (GO) 
annotation analysis was performed to amend informa- 
tion of biological functions [44], including: (1) Enzyme 
Commission number (EC), (2) InterPro scan, (3) 
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annotation augmentation through the Second Layer 
concept (ANNEX) and (4) summarizing annotation 
results by creating a subset of the gene ontology voca- 
bulary encompassing key ontological terms (GOslim). 
Additionally, direct NCBI blasts were performed for 
selected genes of special interest. 
Microarray analysis of the short-term samples (3 and 7 
days) 

Hybridization of the microarray slides was performed 
using the MAUI Hybridization system (MAUI Mixer SC 
Hybridization Chamber mixers, Biomicro Cat. no. 
02A00830) using the ProntolUniversal Hybridization Kit. 
Slides were scanned in an Axon GenePix 4000B dual 
wavelength scanner using the GenePix 6.0 software. 

Initial within-slide normalization (lowess-normaliza- 
tion) was performed by the Acuity 4.0 software (Input 
parameter: Smoothing = 0.19 (chosen by graphical/visual 
evaluation), Iterations = 3 (default), Delta = 0.01 
(default)). For further processing and statistical analysis, 
data was then exported to Excel. Filtering of the data set 
started with excluding low-quality features according to 
the "flagging" performed in Acuity 4.0. Genes were 
defined "absent" and excluded from further analysis, 
when the fluorescence intensity of one or both channels 
was less than 30% higher than the background fluores- 
cence intensity. Also, transcripts with a variance of 
more than 20% in the 4 technical replicates included in 
each slide were not included in the data set. Gene 
expression was calculated as the log2 of the ratio of the 
fluorescence intensity of the C0 2 -treatment cDNA to 
the fluorescence intensity of the control cDNA (log 2 - 
ratio = F635/F532) and as the median and median 
deviation (MD) of each 6 replicate microarray slides per 
block (gill7day3, gill7day7, gill9day3, gill9day7). 

Quantitative polymerase chain reaction (qRT-PCR) 
RNA extraction 

RNA from whole left gills #9 (only day 7, control C0 2 , n 
= 6; elevated C0 2 , n = 6) from the second 1 week 
experiment (2010) and from whole left gills #7 (control 
C0 2 , n = 9; elevated C0 2 , n = 9) of the 11 week experi- 
ment was extracted using the RNeasy Midi-Kit from 
Qiagen (#75144). Gills were first homogenized with an 
OMNI international TH/G7-195STW for 50 s in 4 mL 
buffer RLT from the kit. Extraction then followed the 
manufacturer's description. 
Primers 

A number of transcripts shown to undergo changes in 
the 1 week microarray experiment and/or known to be 
involved in acid-base regulation were selected for analy- 
sis by qRT-PCR, according to their relevant function 
and statistical significance. Primers were designed 
according to the database previously generated by D.W. 
Towle [37], using PrimerExpress software (Version 2.0, 



Applied Biosystems). Primer quality and specificity was 
tested by blasting with the EST library for Carcinus 
maenas from NCBI. Forward and reverse primers were 
generated with the following parameters: GC content 
30-80%, primer length 9-40 bp (opt. 20 bp), Amplicon 
melting temperature 0-85°C, Amplicon length 70-150 bp 
(Additional File 1 Table 2). Performance and efficiency 
of each primer pair was tested in triplets in a cDNA 
dilution series (1:40, 1:80, 1:160, 1:320, 1:640). Perfor- 
mance was evaluated as suitable if the R 2 of the linear 
regression of the dilution series was > 0.98 and no back- 
ground noise (unspecific binding) was detected. A melt- 
ing curve analysis was performed for each reaction to 
ensure a single PCR amplicon. Efficiency was calculated 
as follows: 

E = 10 (1/s) , with S being the slope of the linear 
regression of the dilution series. Only primer pairs with 
an efficiency between 1.9 - 2.1 were chosen for analysis. 
qRT-PCR 

Following DNase digestion of 2 ug total RNA with the 
DNAfree™ Kit from Ambion (#AM1906), 0.4 ug DNA- 
free RNA was then transcribed into cDNA using the 
High Capacity cDNA RT Kit from Applied Biosystems 
(Part no. 4368814). Real-time PCR was performed on 80- 
fold diluted cDNA in 96-well plates using the FastSYBR®" 
Green Master Mix (Part no. 4385612) and the StepOne- 
Plus™ Real-Time PCR System from Applied Biosystems 
(PCR-conditions: 20 s at 95°C, 40 cycles of 5 s at 95°C 
and 30 s at 62°C). On each plate, one control cDNA and 
one experimental cDNA were analyzed with the same 
primer pairs, as well as the respective DNAfree RNA for 
each of the two samples as control for DNA contamina- 
tion, and arginine kinase as housekeeping gene to elimi- 
nate within-plate variation. Arginine kinase was chosen 
as housekeeping gene according to [36] and was tested to 
be appropriate (stability value = 0.007) with NormFinder 
([45], available at http://www.mdl.dk/publicationsnorm- 
finder.htm). Only those samples were included in the 
analysis that showed no or low background in the DNA- 
free RNA control compared to the cDNA samples (differ- 
ence in cycles > 9 » 0.3%). Gene expression was 
calculated based on the C T -threshold. First, absolute 
quantities Q x of all genes were calculated as 

Qx{gene) = (E) Ct , followed by the calculation of nor- 
malized quantities 

QN(gene) = — , with Q x (hk gene) being the 

^>X{HKgene) 

absolute quantity of the housekeeping gene arginine 
kinase. For matters of consistency with the microarray 
data, gene expression resulting from qRT-PCR analysis 
was also calculated as the log 2 -ratio of normalized gene 
quantity in elevated j?C0 2 treated animals to normalized 
gene quantity in control pC0 2 animals. 
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Table 2 Enrichment analysis in the short-term hypercapnia microarray study on Carcinus maenas 



Test-set 


GO 


name 


FDR 


Over-/Under- representation 


down day 3 


bU:UUbU/94 


regulation of cellular process 


0.00 


under 


down day 3 


bU:000/l ob 


signal transduction 


0.00 


under 


down day 3 


bU:Ou2o(JoO 


signal transmission 


0.00 


under 


down day 3 




signaling process 


0.00 


under 


down day 3 




catalytic activity 


0.00 


under 


UUvvl 1 Udy 3 


uU.UUJU/ 


ItrLJ UldUUI I (J! UIUIUL] ILdl piULtrii 


0 04 




up day 3 


no enrichment 








down day 7 


no enrichment 








up day 7 


no enrichment 








up gill 7 


GO:0007010 


cytoskeleton organization 


0.05 


under 


down gill 7 


GO:0005623 


cell 


0.04 


over 


up gill 9 


GO:0005198 


structural molecule activity 


0.03 


over 


down gill 9 


no enrichment 









Results of the enrichment analysis (Fisher's Exact Test with multiple testing correction; FDR < 0.05) as implemented in Blast2GO [43]. Tested were subsets of all 
significantly regulated transcripts as identified by sign test (up- and down-regulated separately for each gill and day or each combination of both). Reference-set 
in each analysis was the complete microarray sequence Blast2GO database generated from the contig list. GO - Gene Ontology term, name = GO name of the 
functional group, FDR = false discovery rate for correction of multiple testing, OverVUnderrepresentation = a GO term is considered over-/under-represented if it 
appears significantly more often/less often in the reference-set than it does in the test-set. 



Statistics 

Because of the high biological variance within replicate 
samples, non-parametric statistical tests were chosen. 
To identify significant changes in gene expression, a 
two-sided sign test with a = 0.05 was performed on the 
fluorescence intensities of the filtered and lowess-nor- 
malized data set for each of the four experimental 
blocks (gill7day3, gill7day7, gill9day3, gill9day7). For 
each individual transcript in each block, the fluorescence 
intensities for the F635-dye, representing the animals 
being exposed to elevated pC0 2 , were tested against the 
control animals (F532-dye). The test only allowed for 
including those transcripts that were represented by 5 
or 6 replicates. 

Sign test was also applied on the experimental vs. con- 
trol normalized gene quantities of the qRT-PCR. 

In order to detect differences between gills and days of 
the short-term exposed animals, Wilcoxon's matched 
pairs test as implemented in STATISTICA 8 was 
applied with a = 0.05. Prior to testing, microarray values 
were transformed into a matrix with 1 = significantly 
up-regulated, -1 = significantly down-regulated, 0 = not 
significantly regulated transcript. Data of each sampling 
day and gill were then tested against each other. 

Additionally, intensity values for both fluorescence 
dyes of the processed microarray data (lowess-normal- 
ized and filtered as described above) was read into the R 
software environment (version 2.13.1, http://www.R-pro- 
ject.org), using the R/maanova package (version 2.22.0, 
[46]). As the package does not allow for missing data, 
only transcripts represented by all 6 replicates could be 



tested by this method. A linear mixed model 

yijk = |x + C02i + Gj + Dk + eijk 

was fitted to the data, with yijk as the normalized- 
transformed gene expression, u as the group mean, 
C02i as the effect of the C0 2 level, Gj as the effect of 
ith gill, Dk as the effect of jth day, and eijk as the sample 
effect (random error). The significance of effects of the 
fixed factors "C0 2 treatment level", "gill" and "day" on 
the differential expression of genes was tested by appro- 
priate contrasts in F-tests between groups in a multifac- 
torial ANOVA design. The F-tests were calculated with 
the James-Stein shrinkage estimate (Fs), incorporating 
shrinkage estimates of variance components [47]. P 
values were calculated by performing 1,000 permuta- 
tions of samples to break their association to expression 
values, then corrected for multiple comparisons by false 
discovery rate transformation (FDR), using the qvalue 
package (jsFDR, version 1.26.0 [48]) with a 5% FDR 
cutoff. 

To identify overall affected pathways, enrichment 
analysis (FDR < 0.05) was performed using Fisher's 
exact test as implemented in Blast2GO [43]. All signifi- 
cantly regulated genes identified by sign test, divided 
into subsets for each gill and day, were used as test- 
sets and tested against a complete microarray sequence 
Blast2GO database generated from the contig list. It 
has to be considered that Carcinus maenas is a non- 
model organism and that therefore the annotation of 
gene functions/Gene ontology (GO) terms is not 
comprehensive. 
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Results 

General findings of the microarray analysis 

Using variance-based analysis utilizing the R/maanova 
package in R, we only identified the factor "C0 2 " to 
have a significant effect on gene expression, while fac- 
tors "Gill" and "Day" did not significantly influence 
expression patterns. 956 out of 3,720 tested transcripts 
were found to be differentially expressed between the 
two C0 2 levels. This accounts for 26% and is compar- 
able to the result of the sign test (24%, see below). 378 
(40%) of the transcripts identified to be differentially 
expressed for the factor C0 2 in the variance-base 
approach were also identified as significantly regulated 
by sign test, including 11 (65%) of the transcripts of spe- 
cial interest as presented in table 3. The sign test identi- 
fied as many as 1,056 (24%) out of the 4,462 transcripts 
on the microarray to be significantly up- or down-regu- 
lated in at least one gill or at one sampling time point 
during the short-term response to elevated seawater 
pC0 2 (p < 0.05; Figure 2; Additional File 1 Table S3). 
Significant up-regulation was observed in 541, down- 
regulation in 502 genes (13 genes showed a mixed regu- 
lation pattern). However, magnitudes of change were 
generally low: expression of only 84 genes of the signifi- 
cantly regulated transcripts (sign test, accounting for 2% 
of the whole microarray and 8% of the significantly 
altered transcripts) were modulated by a log 2 -ratio > | 
0.200 1 (sign test; equivalent to a > 1.3-fold change). F- 
test identified 16 additional transcripts with a log 2 -ratio 
> 1 0.200 1 . 

Enrichment analysis demonstrated that "structural 
molecule activity" (Gene ontology term (GO):0005198) 
was over-represented in all significantly up-regulated 
transcripts of gill 9 (table 3). Additionally, "cell" 
(GO:0005623) was over-represented in all significantly 
down-regulated transcripts of gill 7. No other over- 
representation of any GO-term was identified. However, 
7 GO-terms were under-represented for different sub- 
sets of all significant regulated transcripts in the same 
analysis. 

Regulation of specific groups 

Stress response Out of 137 identified transcripts asso- 
ciated with the cellular stress response (according to 
[38]), only 30 genes were significantly regulated (sign 
test with p < 0.05), albeit less than 1.1-fold (20 up-, 9 
down-regulated transcripts; Additional File 1 Table SI). 
Variance-based linear modeling identified 23 genes to 
be differentially expressed, including 5 that were identi- 
fied also by sign test. Expression of only 5 out of 26 
heat-shock proteins (HSPs) was significantly altered (2 
up-, 3 down-regulated transcripts) in sign test, while 5 
different HSPs were identified by variance-based analy- 
sis. Additionally, stress-associated GO-terms were not 
found to be significantly enriched in the gills of short- 



term exposed animals (Fisher's exact test). The expres- 
sion level of glutathione peroxidase [GenBank: 
DN551450.1] was also investigated in the long-term 
exposed animals (qRT-PCR analysis) and did not show a 
significant change. 

Structural modification In total, 77 transcripts were 
associated with GO: 0005198 (structural molecule 
activity). Significant up-regulation of 13 genes (sign 
test, 10 of them also identified by F-test) from this 
group resulted in the over-representation of this GO- 
term in gill 9 (Additional File 1 Table S5). These up- 
regulated genes are eight ribosomal proteins, a cuticu- 
lar protein, two alpha-tubulins, a keratin-associated 
protein and a vitelline membrane outer layer protein. 
The identification of a 1.4-fold up-regulated transcript, 
encoding for a multispanning endomembrane protein 
of the transmembrane 9 superfamily (protein member 
4 = TM9SF4, [GenBank:DN202592.1]), also suggests 
structural modification of gill membranes. It was up- 
regulated in both gills on both sampling days (sign test 
with p < 0.05). This transcript was also identified to be 
up-regulated in the low salinity study [37]. Addition- 
ally, two other transmembrane proteins were identified 
to be strongly down-regulated (sign test with p < 0.05): 
a tetraspanin-like 8 protein ([GenBank:DV944030.1], 
included in the GO-class "cell") was found to be 
strongly down-regulated during short-term hypercap- 
nia (1.3-fold) on day 3 and slightly during low salinity 
exposure. A putative Integrin-alpha-7 [GenBank: 
DY307809.1] was identified to be strongly down-regu- 
lated (1.3-fold, sign test with p < 0.05) on day 7 (table 
3, Additional File 1 Table S5). 

Specialization of the gills Although no effect of the 
factor "gill" on individual gene expression was identified 
by variance-based linear modeling, Wilcoxon's matched 
pairs test identified a significant difference in the overall 
gene expression patterns in the different experimental 
blocks (gill7day3, gill7day7, gill9day3, gill9day7; p < 0.05; 
Figure 2, Figure 3A), based on total numbers of up-, 
down-, or not regulated transcripts. Fewer genes were 
significantly regulated in gill 9 after both 3 and 7 days 
of exposure than in gill 7 (sign test). In total, 756 genes 
were significantly regulated in gill 7, whereas only 510 
in gill 9. Upon hypercapnia stress, gill 7 in general 
shows stronger responses to hypercapnia than gill 9. For 
example, the strongest overall significant change in gene 
expression was observed in gill 7 of a senescence-asso- 
ciated protein (sign test and F-test with p < 0.05, [Gen- 
Bank:DV944570.1]), a gene that was not altered in gill 9, 
regarding both, day 3 and 7. On the other hand, gene 
expression of a calcium-activated chloride channel [Gen- 
Bank:DW584691.1] was only significantly and strongly 
affected in gill 9 on both days of the short-term expo- 
sure (sign test and F-test with p < 0.05). Additionally, 



Table 3 Regulation of transcripts of special interest in the short-term hypercapnia study on Carcinus maenas 
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Included are the 14 most strongly regulated transcripts (as identified by sign test and F-test), as well as several transcripts which have a potential function during structural modifications or in ion- and acid/base 
regulation. Transcripts are ordered according to the rank in absolute regulation. Accession numbers (ACC. no.) refer to the ESTs generated for the C. maenas microarray by Towle et al. [37] in the database GenBank 
(NCBI). Rank 1 = highest regulation; bold and underlined = significant regulation relative to the control identified by sign test; bold = identified to be differentially expressed by both tests; italic = identified to be 
differentially expressed only by F-test; "-" = a rank was not assigned in case the regulation was less than log 2 -ratio = |0.20|. 
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Figure 2 Differentiation of the sum of all significantly 
regulated transcripts for each gill and day of Carcinus maenas 
upon short-term exposure to hypercapnia For each of the 
experimental blocks of the short-term hypercapnia experiment 
(gill7day3, gill7day7, gill9day3 and gill9day7), up- and down- 
regulated transcripts were counted and represented as bars. 
Additionally, a matrix was generated for all significantly regulated 
transcripts (sign test) for each block with 1 = significantly up- 
regulated, 0 = not significantly regulated, -1 = significantly down- 
regulated. Statistical comparison of the four experimental blocks 
was done applying Wilcoxon's matched pairs test on the matrices. 
No significant difference was found in gill 7 comparing the two 
sampling points, while gill 9 on day 3 differed significantly 
compared to gill 9 on day 7 (p < 0.05). Additionally, gill 9 in genera 
differed significantly in comparison to gill 7 on both days (p < 0.05). 
Different letters a, b, c denote these significant difference. 



according to the enrichment analysis, structural changes 
in gill epithelia are linked mainly to gill 9. 
Comparison with gene expression patterns in response 
to low salinity When compared to gene expression 
changes observed after 7 days of acclimation to a 2-fold 
reduction in salinity (from S = 32 to S = 15; [37]), 
short- and long-term acidification (10-fold increase in 
j?C0 2 ) resulted in smaller changes in expression levels 
in fewer transcripts: only 8 transcripts were regulated > 
1.5-fold (log 2 -ratio > |0.60|) in the present study, while 
533 transcripts were regulated > 1.5-fold in the low sali- 
nity acclimation study; maximum expression changes of 
up to 14-fold (log 2 -ratio = 3.8) were observed in the low 
salinity acclimation study (Figure 3B, 4). Nevertheless, a 
considerable amount of transcripts were regulated dur- 
ing acclimation to both abiotic stressors (table 2). 
Comparison of qRT-PCR and microarray analysis In 
the long-term experiment, no significant changes in 
expression were found, likely due to the high biological 
variance within samples (coefficients of variation of the 
ratio 400 Pa/39 Pa CO a between 15 to 79%, Figure 5). 
In the short-term experiment, 7 out of the 8 tested 
genes responded in the same direction as in the micro- 
array analysis (Additional File 2 Figure SI). In one case 




significant in gill 7, day 3 
significant in gill 9, day 3 
significant in gill 7, day 7 
significant in gill 9, day 7 
not significant 



CD 
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log2-ratio [hypercapnia; gill 7, both days] 




gj -3 -2 -1 0 1 2 3 

o 

log2-ratio [salinity; gill 8, day 7] 

Figure 3 Comparison of effects in transcript regulation in gills 
of Carcinus maenas upon short-term exposure to hypercapnia 
and low salinity. (A) Regression analysis of gill 9 vs. gill 7 (both 
days) in the short-term hypercapnia experiment and (B) of gill 7/day 
7 (short-term hypercapnia experiment, this study) vs. gill 8/day 7 
(salinity dilution experiment by Towle et al. [37]). The long-dashed 
black line represents a theoretical regression ratio of 1:1 (equal 
influence of both factors), while the black line shows the actual 
regression. Transcripts outside the short-dashed horizontal and 
vertical lines are more than 2-fold regulated. Please note the 
different scaling in (A) and (B). Numbers are according to the gene's 
rank in table 3, indicating transcripts discussed in further detail in 
the manuscript. (1) senescence-associated protein [GenBank: 
DV944570.1], (3) putative Syntaxin binding protein 2 [GenBank: 
DN161290.1], (4) 23S ribosomal RNA gene (a) [GenBank:DN796131.1], 
(5) calcium-activated chloride channel [GenBank:DW584691.1], (6) 
23S ribosomal RNA gene (b) [GenBank:DN635040.1], (13) 
multispanning membrane protein (transmembrane 9 superfamily 
protein member 4 [GenBank:DN202592.1], (25) potential integrin 
alpha-7 [GenBank:DY307809.1]. 



(glycosyl-phosphatidylinositol-linked carbonic anhydrase 
VII [GenBank:DN739347.1]), the transcript was slightly, 
but significantly down-regulated in the microarray ana- 
lysis, whereas it was up-regulated in the qRT-PCR ana- 
lysis. Variation in both, qRT-PCR and in the microarray 
experiment was high. 
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Figure 4 Comparison of regulation frequencies in Carcinus 
maenas gills in salinity and C0 2 experiments. (A) Frequency 
distribution of all gene expression levels in gill 7 and 9 on day 7 as 
response to short-term hypercapnia (this study). (B) Frequency 
distribution of all gene expression levels in gill 8 on day 7 as 
response to low salinity (study of Towle et al. [37]). Please note the 
different scaling of the y-axis in (A) and (B). Exposure to low salinity 
(B) clearly resulted in more genes to be changed on a higher level 
compared to hypercapnia (A), where levels of change were mainly 
between a log 2 -ratio of -0.1 25 - 0.1 25. (C) Numbers of genes found 
to be regulated in expression simultaneously in response to a 
salinity dilution, as well as short-term hypercapnia (only transcripts 
with a log 2 -ratio > |0.200| were considered). 
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Figure 5 Change of expression levels of distinct transcripts in 
gill 7 as response of Carcinus maenas to long-term exposure 
(11 weeks) to hypercapnia analysed by quantitative real-time 
polymerase chain reaction. Bars represent the median log 2 -ratio 
and median deviation (error bars) of the respective gene. Numbers 
and abbreviations refer to primers and genes as used in the qRT- 
PCR (see Additional File 1 Table S2 for primer sequences and full 
names). No statistical significance was detected (sign-test with p < 
0.05), but several transcripts exhibit a tendency to be up- or down- 
regulated, e.g. chloride associated transporters (no. 16, 18, 21). 



Most strongly affected transcripts The most down- 
regulated gene coded for a senescence-associated pro- 
tein [GenBank:DV944570.1] in gill 7 on day 7 with a 
median log 2 -ratio of -1.16 ± 0.64 (2.2-fold change, sign 
test and F-test with p < 0.05, table 3). Its expression was 
not altered in the low-salinity experiment. An unknown 
transcript [GenBank:DN161290.1], suggested to encode 
for Syntaxin binding protein 2 by Towle et al. [37]), was 
the most strongly up-regulated transcript in gill 7 on 
day 7 with a median log 2 -ratio of +0.88 ± 0.36 (1.9-fold 
change, sign test and F-test with p < 0.05). This tran- 
script was also the most strongly up-regulated gene in 
gill 8 in response to low salinity transfer [37]. Pgl pro- 
tein [GenBank:DW250369.1], only identified by F-test, 
was another highly up-regulated transcript on day 3, 
while it showed strong down-regulation on day 7. In gill 
9, the most strongly down-regulated gene was identified 
on day 7 as a 23S ribosomal RNA gene [GenBank: 
DN796131.1] with a median log 2 -ratio of -0.87 ± 0.46 
(1.9-fold change, sign test and F-test with p < 0.05). 
This gene was strongly up-regulated during the low-sali- 
nity experiment. Additionally, a different 23S ribosomal 
gene was found to be equally strongly affected on day 7 
in gill 7 ([GenBank:DN635040.1], log 2 -ratio = 0.80 ± 
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0.22, 1.7-fold change, sign test with p < 0.05, not 
included in F-test). The most strongly up-regulated 
transcript in gill 9 at day 3 was a calcium-activated 
chloride channel [GenBank:DW584691.1] with a median 
log 2 -ratio of +0.81 ± 0.34 (1.7-fold change, sign test and 
F-test with p < 0.05). It was also the most highly up- 
regulated transcript in the long-term exposed animals 
(log 2 -ratio = 1.12 ± 0.88, 2.2-fold change), although not 
tested significant due to high variance between repli- 
cates. Additionally, rRNA intron-encoded homing endo- 
nuclease [GenBank:DV944193.1] was identified by F-test 
as being differentially expressed. Besides 4 unknown, 
but strongly changed transcripts (rank 8-11, table 2), a 
thymosin isoform was identified to be strongly down- 
regulated in gill 9 on day 3 (log 2 -ratio = -0.46 ± 0.14, 
1.4-fold change, sign test with p < 0.05). In addition, a 
hyperpolarization-activated cyclic nucleotide-gated 
potassium channel (HCN2, [GenBank: DW250260.1]), 
was significantly down-regulated in both gills on day 7 
(gill 7 with a log 2 -ratio = -0.46 ± 0.24/1.4-fold, and gill 
9 with -0.14 ± 0.08/1.1-fold, sign test and F-test with p 
< 0.05) with a tendency for down-regulation in both 
gills on day 3. In the low-salinity study, it was only 
slightly down-regulated [37]. 

Acid/base and ion regulation None of the acid-base 
regulatory candidate genes were regulated more than 
1.3-fold in our microarray experiment at days 3 and 7. 
Of 9 candidate genes tested on the microarray, only 2 
were significantly regulated with respect to controls. 
Among them were a CI7HCO3" exchanger ([GenBank: 
CX994129.1], -1.1-fold, sign test and F-test with p < 
0.05) and the glycosyl-phosphatidylinositol-linked carbo- 
nic anhydrase VII ([GenBank:DN739347.1], -1.1-fold, 
sign test and F-test with p < 0.05). In contrast, during 
acclimation to reduced salinity, Towle et al [37] found 5 
of this 9 acid-base regulatory candidate genes to be 
altered by more than 1.3-fold, including the cytoplasmic 
carbonic anhydrase ([GenBank:DW251095.1], + 3.4- 
fold), the Na + /K + -ATPase alpha subunit ([GenBank: 
DW250552.1], + 2.6-fold), a member of the anion/bicar- 
bonate transporter family (abts-3) ([GenBank: 
DN202373.1], similarity to a Na7HC0 3 ~ transporter 
(SLC4A11), - 2.4-fold), the glycosyl-phosphatidylinositol- 
linked carbonic anhydrase VII ([GenBank:DN739347.1], 
- 1.9-fold) and the C17HCGY exchanger (SLC4A1; 
[GenBank:CX994129.1], + 1.3-fold). 

In the long- term acclimation study (11 weeks), no 
significant changes in expression could be found, but 
variances were extremely high (coefficients of variation 
of the ratio 400 Pa/39 Pa CO a ranged from 15 to 79%). 
Only in the case of the sodium-hydrogen exchanger 
NHE3, variance was smaller (1.3%). 



Discussion 

Green crabs, especially those of the Baltic Sea, are 
known to be effective ion and acid-base regulators 
[25,49,50]. When exposed to hypercapnia, Carcinus 
maenas is characterized by a rapid pH e compensatory 
response that involves the accumulation of bicarbonate 
in its hemolymph ([17]). This high HC0 3 ~ level is sus- 
tained even over a period of 11 weeks (Appelhans et al., 
under review). The enrichment of HCO3' in the crabs' 
hemolymph may be mediated by altered ion regulation 
processes in the gills and thus can be expected to also 
leave a footprint in the gills' transcriptome. However, 
long-term hypercapnia incubations may additionally 
result in other physiological responses, pointing to 
trade-offs in energy allocation in favor of extra-/intracel- 
lular ion homeostasis vs. protein metabolism and 
growth, as has been concluded for molluscs, teleost fish 
and echinoderms [13,23,51-53]. We hypothesized that 
an increased demand for ion transport and modifica- 
tions in epithelial C0 2 permeability would influence 
expression patterns of the respective gene transcripts in 
gills in both, short- and long-term hypercapnia 
scenarios. 

Variance-based analysis of the microarray identified a 
significant effect of pC0 2 (53 vs. 440 Pa) on 26% of the 
transcripts. However, expression profiles in this study 
revealed that short-term hypercapnia (10-fold increase 
in seawater j?C0 2 ) does not act as a strong stressor on 
gill tissues. Expression changes were moderate and it 
has to be considered that some of the observed changes 
are random noise, given also that only few significant 
effects persisted over both time points (3 and 7 days) 
and considering the partially differing results from both 
statistical tests applied (40% of trancripts identified to 
be differentially expressed by both, sign test and var- 
iance-based analysis). Only relatively few transcripts (2% 
of all tested transcripts) were differentially regulated 
more than 1.3-fold, with maximum observed changes of 
2.2-fold. Even these changes are much lower than those 
elicited by a 2-fold decrease in salinity (transfer from S 
= 32 to S = 15), where maximum observed changes in 
expression were 14-fold [37]. A negative effect of hyper- 
capnia on transcripts related to energy metabolism, as 
has been shown in sea urchin larvae [53], was not 
observed in our study. Transcripts known to be involved 
in acid-base responses were not regulated strongly in 
response to elevated pC0 2 - Instead, new acid-base regu- 
latory candidate genes were identified to be differentially 
expressed, including a calcium activated chloride chan- 
nel, a multispanning membrane protein (transmembrane 
9 superfamily protein member 4), a hyperpolarization 
activated cyclic nucleotide-gated potassium channel-2, 
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tetraspanin 3 (transmembrane super family 4 member 8) 
and a potential integrin alpha-7. 

Short-term hypercapnia does not lead to a pronounced 
stress response 

Only few stress-associated genes were differentially 
regulated in in gills in response to short-term hypercap- 
nia. Therefore, short-term hypercapnia does not evoke a 
pronounced cellular stress-response (CSR) in the gills of 
adult Carcinus maenas from the Baltic Sea. A lack of 
pronounced changes in the expression of heat-shock 
proteins supports this interpretation, as does the finding 
that in general only small expression level changes 
occured. Heat-shock proteins (HSPs) are highly con- 
served among species and involved in acute stress 
responses towards various stressors [54]. Usually, cellu- 
lar responses through HSP occur rapidly and transiently 
[55,56]. Thus, it cannot be excluded that a CSR 
occurred and ended before our first sampling after 3 
days. However, the observation that a hypoosmotic 
acclimation with its much stronger changes in gene 
expression patterns also did not elicit a pronounced 
CSR (even within the first hours of exposure; [37]) cor- 
roborates the interpretation that the applied changes in 
seawater pC0 2 do not trigger a prolonged or pro- 
nounced CSR in gills of C. maenas. As focus has been 
laid upon the gills, it should be considered that hyper- 
capnia may indeed provoke a CSR in other organs. 

Although magnitudes of change in gene expression 
were comparatively low (e.g. compared to [37,53]), we 
found indications for a re-arrangement of the green 
crabs' gill epithelia and support for specialization of 
each single posterior gill. 

Re-arrangement of the gill epithelium and epithelial cell 
membranes 

Enrichment analysis supports that a reorganization of 
gill epithelia occurs during the initial stages (< 7 days) 
of hypercapnia acclimation. The GO-term "structural 
modification" is the only biological process that was 
enriched within up-regulated genes in response to 
hypercapnia. During salinity acclimation, posterior gills 
of C. maenas also undergo structural modification 
through extension of the apical plasma membrane 
infolding system of the thick prismatic salt-transporting 
epithelium and through an increase of the subcuticular 
space [50]. Although no genes encoding for cellular 
junction proteins such as claudins, occludins, cadherins 
or selectins were tested on the microarray used in the 
present study, our analysis identified several other genes 
potentially involved in structural rearrangements. A 
member of the tetraspanin family was found to be 
strongly down-regulated during short-term hypercapnia 
and slightly during low salinity acclimation [37]. 



Tetraspanins are a group of four-transmembrane- 
domain proteins that are expressed in epithelia and 
known to be involved in diverse cellular processes (e.g. 
morphogenetic re-organization of monolayers of epithe- 
lial cells [57]). They are generally described as molecular 
'facilitators' or 'organizers' at the plasma membrane 
[58,59]. Furthermore, integrin-tetraspanin complexes 
play an important role in cell-cell-adhesion at cellular 
junctions [56,59,60]. A putative Integrin-alpha-7 was 
identified to be significantly and strongly down-regu- 
lated in the short-term acidification experiment. Due to 
their essential role in cell adhesion and cell-cell commu- 
nication, biochemical functions of integrins are likely to 
be highly conserved in metazoans [61]. The common 
and distinct regulation of a tetraspanin and an integrin 
indicate that the possible complex of both could play a 
role during re-arrangement of the gill. Another strongly 
up-regulated transcript encodes for a multispanning 
endomembrane protein of the transmembrane 9 super- 
family (protein member 4 = TM9SF4), also known as 
p76 in humans. TM9SF4 plays an important role in cel- 
lular adhesion, membrane reconfiguration and vesicle 
mediated transport [62,63]. It thus can be hypothesized 
that exposure to elevated seawater pC0 2 has an influ- 
ence on the membrane composition of gill epithelial 
cells, and on the cell composition of the epithelium of 
the gills of C. maenas. This warrants detailed studies on 
structural changes of the gills in response to hypercap- 
nia (eg. light or electron microscopy and (immune-) his- 
tochemical investigations), including studies on the 
involvement of the three membrane proteins mentioned 
above. 

Specialization of the gills 

As described above, only the three posterior gills 7 to 9 
were found to be involved in ion- and acid-base regula- 
tion in C. maenas. The anterior six gill pairs are primar- 
ily important for gas exchange [17,26,27]. However, a 
more individualized specialization of each of the poster- 
ior gills has been shown for other species with respect 
to salinity changes [64,65]. For C. maenas, Siebers et al. 
[26] detected a salinity-dependent activation of Na + /K 
+ -ATPase with increasing activity at decreasing salinities 
mainly in the posterior gills, but only subtle differences 
between the individual gills 7-9. Henry et al. [66] made 
a corresponding observation with respect to carbonic 
anhydrase. However, based on the distribution of V-type 
H + -ATPase in the posterior gills of several (intertidal) 
crab species, Tsai and Lin [65] showed that the func- 
tional differentiation in crab gills generally is not only 
between anterior and posterior, but also within indivi- 
dual gill lamellae. While the variance-based analysis of 
the data set suggested no significant effect of the factor 
"gill" on gene expression levels, enrichment analysis 
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indicated different and specific responses to hypercapnia 
acclimation between posterior gills 7 and 9. Gill 7 
responded stronger to hypercapnia than gill 9 with 
respect to both, the number of genes affected and the 
overall magnitude of change in gene expression. Most 
transcripts were either significantly regulated for one or 
the other gill, only 11% were regulated in parallel. The 
result of the enrichment analysis further suggested that 
structural changes are mainly associated with gill 9. 
However, the conflicting results of the two statistical 
methods used demand more detailed investigations on 
the cell ultrastructural level to substantiate potential dif- 
ferences in function between gills. 

Comparison with gene expression levels in response 
towards low salinity 

When responses to hypercapnia (this study) and hypoos- 
motic acclimation [37] are compared, it becomes 
obvious that hyposmotic acclimation to a 2-fold reduc- 
tion in salinity results in far larger expression changes 
than acclimation to a 10-fold increase in seawater pC0 2 . 
Nevertheless, acute hyposmotic acclimation also does 
not lead to a pronounced CSR [37]. C. maenas from the 
western Baltic Sea show an increased ability for hyperos- 
moregulation when exposed to low salinities (salinity < 
20) than animals from the more saline North Sea (sali- 
nity > 30, [47]). Consequently, they also must possess a 
very high ion regulatory capacity and as salinity was low 
in our experiment (S = 15), the gill ion regulatory 
machinery was working at a comparatively high load. 
Microarray analysis only allows detection of relative 
changes in gene expression. If a certain transcript is 
already highly expressed, even a small relative change 
on the level of transcript expression could result in a 
highly effective regulatory capacity on the protein level. 
Even small changes may allow for successful acclima- 
tion. Furthermore, the strong fluctuations (both in rate 
and magnitude) in salinity, pC0 2 and temperature 
observed in Kiel Fjord [19], might have led to an adapta- 
tion towards excess ion regulatory capacity that can be 
recruited upon demand. This recruitment might take 
place on the post-transcriptional level [23,67] and would 
therefore remain undetected in our microarray analysis. 

Old and new candidate genes for hypercapnia 
acclimation in crustaceans 

In addition to the above discussed potential role of the 
multispanning membrane protein TM9SF4 in structural 
re-arrangement of gill epithelia, its hypercapnia induced 
up-regulation might be relevant for cellular acid-base 
regulation. TM9SF4 is known to participate in vesicular 
transport [68] and Schimmoller et al. [69] suggested its 
association with endosomes (acidic compartments/vesi- 
cles in mammalian cells). In posterior C. maenas gills, 



intracellular vesicles were postulated to be involved in 
cellular acid-base regulation via V-H + -ATPases [70]. 

A calcium-activated chloride channel (CaCC) was 
strongly up-regulated in gill 9 on day 3 and it was also 
up-regulated after 11 weeks of hypercapnia. Beside 
others, CaCCs have been shown to play a key role in 
epithelial secretion [71,72], but a high variability within 
this class of channels with respect to physiological 
roles and mechanisms of regulation was observed 
[73,74]. A member of a CaCC subfamily has also been 
shown to act in cell-cell adhesion through interaction 
with an integrin (e.g. in lung cancer [67], reviewed 
[75]). Chloride channels have been hypothesized to be 
situated in the basolateral membrane of epithelial cells 
in gills of osmoregulating crabs [32,33]. The identified 
CaCC as described above might be an important can- 
didate gene in CI" regulation. In order to achieve elec- 
troneutrality, CI" typically is the counter-ion of HC0 3 " 
during the extracellular pH regulatory reaction. Extra- 
cellular HCGV accumulation is probably enabled by 
C17HCGV exchangers [20]. While respective C17 
HCGV exchangers (AE) have been postulated to be 
situated in the apical membrane in crustaceans, mole- 
cular identification and/or biochemical characterization 
is still lacking. On the other hand, a basolateral-situ- 
ated AE can be discussed to play a role in pH and 
volume regulation [33]. In acid secretion, a respective 
exchanger is postulated to transport HC0 3 " ions from 
the cell into the hemolymph in exchange for CI" ions 
and therefore is argued to sit in the basolateral mem- 
brane. As has been postulated in fish, other transpor- 
ters in close proximity to a C17HC0 3 " exchanger can 
be discussed to favor the electroneutral exchange of 
CI" and HC0 3 " against an unfavorable CI" gradient 
[76]. A Cl"-channel like the identified CaCC could be 
an important additional player and facilitator in this 
CI /HCO3" exchange. 

Additionally, a hyperpolarization-activated cyclic 
nucleotide-gated potassium channel (HCN) was signifi- 
cantly down-regulated in the short-term hypercapnia 
study. So far, HCNs (1+4) have been shown to be situ- 
ated in the plasma membrane of vallate papillea taste 
cells in rat tongue. Those transporters/receptors are 
associated with the basolateral membrane and play a 
role in response to sour stimuli (extracellular protons) 
by mediating an inwardly directed current [77]. On the 
other hand, lowered intracellular pH has been shown to 
lead to a decreased opening speed of the channel in tha- 
lamocortical neurons of the rat ventrobasal thalamic 
complex [78]. Thus, an altered extracellular acid-base 
status might interact with the function of this protein. 

Future studies should characterize these candidate 
genes with respect to localization and function. In the 
case of the CaCC and HCN, electrophysiological 
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experiments could reveal in which way these transpor- 
ters mediate ion fluxes across the gill epithelium. 

Several transporters and channels shown to be 
involved in ion or acid-base regulation during hypercap- 
nia acclimation in other studies on diverse marine 
organisms [33,51,79] were not affected in C. maenas gill 
tissue. The sodium pump, Na + /K + -ATPase (NKA), is 
crucial for the maintenance of ion gradients that drive 
acid-base regulation [80]. This primary active transpor- 
ter is a key player in establishing the characteristic ion 
gradient that is used by many secondary active transpor- 
ters. NKA transcript and protein level, as well as activity 
increased in response to decreased salinities in diverse 
crabs, including C. maenas [C. maenas: 37, 69; others: 
64, 65, 81]. During hypercapnia acclimation, Deigweiher 
et al. [52] documented that in teleost fish, NKA mRNA 
concentration decreased initially (day 4), only to increase 
2-fold after 6 weeks of exposure to a pC0 2 of 1 kPa. In 
contrast, O'Donnell et al. [53] found that in sea urchin 
larvae, NKA mRNA was down-regulated in response to 
hypercapnia. Although NKA expression decreased signif- 
icantly in cephalopod embryos and hatchlings exposed 
to a pC0 2 of 0.4 kPa, no change in expression was 
detected in juveniles under comparable conditions [23]. 
We did not find changes in NKA expression in C. mae- 
nas gill tissue in the present study. Only one C17 
HC0 3 '-anion exchanger from the SLC family 4 (mem- 
ber 1; [GenBank:CX994129.1]) was significantly down- 
regulated in the short-term experiment (1.1-fold), while 
a different anion-bicarbonate exchanger, similar to the 
SLC family 4, member 11 [GenBank:DN202373.1]), and 
a vacuolar H + -ATPase [GenBank:DY656042.1] were not 
affected. Another important transporter for acid-base 
regulation, a Na + /H + exchanger (NHE3) was not signifi- 
cantly regulated in C. maenas gill tissue in the long- 
term experiment. C. maenas also possesses two bran- 
chial isoforms of carbonic anhydrase (CA) in its poster- 
ior gills, a membrane-associated and a cytoplasmic form 
[82]. The enzyme catalyzes the highly energy-demanding 
transition of H 2 0 and CO2 to HCO3" and H + and vice 
versa. It is of great importance in osmoregulation (espe- 
cially the cytoplasmic pool [66]), acid-base regulation 
and C0 2 excretion in the gills of crustaceans [71,83]. 
However in this study, no response of CAs was observed 
in short-term hypercapnia experiments, except for a 
slight significant down-regulation of glycosyl-phosphati- 
dylinositol-linked carbonic anhydrase VII [GenBank: 
DN739347.1] in both, gill 7 and 9, on day 7. In agree- 
ment with the results of the short-term experiment, we 
also found no significant expression changes following 
long-term exposure to hypercapnia. 

We suggest that the ion regulatory apparatus of Baltic 
C. maenas already works at a high load due to the 
demands of a hyposmotic habitat with large fluctuations 



in pC0 2 . Therefore, only moderate mRNA expression 
level changes might be necessary to compensate hyper- 
capnia in C. maenas. Compensation of the ion regula- 
tory apparatus might additionally take place on the 
post-transcriptional level or is facilitated by transcripts 
not included on this microarray. It has to be considered 
though, that to some extend, effects might be hidden 
under the experimental error. Nevertheless, those few 
genes for which we have identified relatively strong 
changes in expression levels are particularly interesting. 
They are likely key players of hypercapnia acclimation 
of crustacean gill tissues. 

Conclusions 

The response of Carcinus maenas to elevated seawater 
pC0 2 based on the gill transcriptome does not suggest 
that seawater acidification acts as a strong stressor for 
the western Baltic population of this species. Following 
short-term hypercapnia, the low but specific responses 
of the gills indicate that (1) the response to hypercapnia 
is partially similar to the response to hypoosmotic con- 
ditions; (2) a multispanning membrane protein 
(TM9SF4), a calcium-activated chloride channel (CaCC) 
and a hyperpolarization-activated cyclic nucleotide-gated 
potassium channel (HCN) are strongly regulated in 
response to hypercapnia and may be involved in acid- 
base regulation; (3) posterior gills might respond differ- 
ently to hypercapnia; and (4) a re-configuration of the 
epithelial gill membrane, including the involvement of a 
tetraspanin-integrin complex, might occur. 

Due to the low salinity environment in the Baltic Sea 
in this hypercapnia study, effects on distinct gene 
expression levels to one of the respective abiotic factors 
might be influenced by the other. It is therefore essential 
to study both abiotic factors separately and in more 
detail. Here, potential pathological effects of the long- 
term consequences of elevated blood bicarbonate con- 
centrations on various tissues and organs should also 
deserve particular attention. We also expect that stron- 
ger effects of hypercapnia are to be found on the pro- 
teome and metabolome level. 

Additional material 



Additional File 1: Tables S1-S6. Table SI. Annotation details on 
transcripts for the Carcinus maenas the microarray assay applied in 
this study. Annotation results on transcripts for the Carcinus maenas 
microarray assay applied in this study from a newly performed 
bioinformatic analysis conducted with the internet portal Blast2G0 [43] 
in 2009 by the authors of this study, compared to a first annotation 
provided by Towle et al. [37]. Included are the sequence length of the 
aligned sequence (= seq. length), the number of hits for the transcript (= 
#hits), the mean similarity of the alignment (= mean sim), the number of 
Gene Ontology (GO) terms assigned to that blast hit (= #GOs), a 
description of the assigned GOs (= GOs), Enzyme Codes (= EC), results of 
the InterProScan, the exact hit description of the best hit, as well as the 
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accession number (= hit ACC, E-value, similarity, score, alignment length 
and number of positive matches bases of the best hit. Table S2. Details 
on primers for quantitative real-time polymerase chain reaction 
(qRT-PCR) to assess responses of Carcinus maenas to hypercapnia 

Primer sequences (5' — > 3') and descriptions of the targeted genes used 
in the real-time polymerase chain reaction (qRT-PCR) in the short- and 
long-term hypercapnia experiments on Carcinus maenas. Numbers ('no.') 
are according to the numbers used in Figure 4 and Additional File 2 
Figure SI. Accession numbers (ACC. no.) refer to the ESTs generated for 
the C. maenas microarray by Towle et al. [37] and the database GenBank 
(NCBI). R 2 and efficiency were tested in a qRT-PCR dilution series (for 
details, see Material & Methods). Table S3. Transcripts significantly 
regulated in response to hypercapnia in Carcinus maenas Transcripts 
significantly regulated in response to hypercapnia in Carcinus maenas as 
identified by variance-based linear modelling (F-test) and sign test. Out 
of 1634 genes, 678 were identified by sign test, 578 genes were 
identified by F-test, and 378 were identified to be affected by both 
statistical tests. Transcripts are sorted alphabetically after the accession 
number (ACC. no., database GenBank (NCBI)). Bold and underlined = 
significantly regulated as identified by sign test, p < 0.05; bold = 
significantly regulated as identified by both statistical tests (sign test and 
F test, p < 0.05); italic = differentially expressed as identified by F-Test. 
Accession numbers refer to the ESTs generated for the C. maenas 
microarray by Towle et al. [37]. Values are given as median and median 
deviation (MD) Table S4. Significantly regulated transcripts 
associated with a cellular stress-response in Carcinus maenas when 
exposed to hypercapnia. Significantly regulated transcripts in the short- 
term hypercapnia study on Carcinus maenas associated with a cellular 
stress-response acoording to [38]. Bold and underlined (values) = 
significantly regulated (sign test, p < 0.05), bold (ACC. no.) = transcripts 
identified to be significantly regulated by both tests (sign test and F-test, 
p < 0.05). Accession numbers (ACC. no.) refer to the GenBank database 
(NCBI).Table S5. Significantly up-regulated transcripts in the short- 
term hypercapnia study associated with structural molecule activity 
in gills of Carcinus maenas Transcripts identified to be responsible for 
the over-representation of the GO-term "structural modification" 
(GO:0005198) in Carcinus maenas gills exposed to hypercapnia 
(enrichment analysis with Fisher's exact test). Bold and underlined 
(values) = significantly regulated as identified by sign test (p < 0.05); bold 
(ACC. no.) = identified to be significantly regulatd by both tests (sign test 
and F-test, p < 0.05). Accession numbers (ACCs) refer to the database 
GenBank (NCBI). Table S6. Details on transcripts of special interest 
identified by microarray analysis on gills of Carcinus maenas after 
exposure to short-term hypercapnia. Identification and details on 
transcripts of special interest identified by microarray analysis, after an 
additional NCBI blastx. * = NCBI blastn; ** = see Towle et al. [37] for 
details. Accession numbers (ACC. no.) refer to the database GenBank 
(NCBI). 

Additional File 2: Figure SI. Comparison of distinct transcripts of 
the microarray analysis vs. qRT-PCR for the response of Carcinus 
maenas to short-term hypercapnia. Comparison of the regulation of 
distinct transcripts of gill 9 for the Carcinus maenas response to short- 
term hypercapnia (1 week, April 2009) in the microarray analysis with 
results of the qRT-PCR experiment performed on the respective genes 
from the short-term incubation conducted in April 2010. In 7 of 8 cases, 
both techniques show the same tendency in regulation. Values represent 
median log 2 -ratios with median deviation (error bars). Transcript numbers 
according to Additional File 1 Table S2. Senesc. ass. prot = senescence- 
associated protein, put. syntaxin = putative Syntaxin binding protein 2, 
gCA = glycosyl-phosphatidylinositol-linked carbonic anhydrase VII, prot. 
inh. = hemozyte kazal-type proteinase inhibitor, K-channel = 
hyperpolarization activated cyclic nucleotide-gated potassium channel 2, 
NKA = Na + /K + -ATPase alpha subunit, Cl-channel (Ca-act.) = calcium 
acitvated chloride channel. 
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